Electron Doping of Cuprates via Interfaces with Manganites 
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The electron doping of undoped high- Tc cuprates via the transfer of charge from manganites (or 
other oxides) using heterostructure geometries is here theoretically discussed. This possibility is 
mainly addressed via a detailed analysis of photoemission and diffusion voltage experiments, which 
locate the Fermi level of manganites above the bottom of the upper Hubbard band of some cuprate 
parent compounds. A diagram with the relative location of Fermi levels and gaps for several oxides 
is presented. The procedure discussed here is generic, allowing for the qualitative prediction of the 
charge flow direction at several oxide interfaces. The addition of electrons to antiferromagnetic Cu 
oxides may lead to a superconducting state at the interface with minimal quenched disorder. Model 
calculations using static and dynamical mean-field theory, supplemented by a Poisson equation 
formalism to address charge redistribution at the interface, support this view. The magnetic state 
of the manganites could be antiferromagnetic or ferromagnetic. The former is better to induce 
superconductivity than the latter, since the spin-polarized charge transfer will be detrimental to 
singlet superconductivity. It is concluded that in spite of the robust Hubbard gaps, the electron 
doping of undoped cuprates at interfaces appears possible, and its realization may open an exciting 
area of research in oxide heterostructures. 

PACS numbers: 73.20.-r, 74.78. Fk, 73.40.-c 



INTRODUCTION 



its study is only now starting [1, 0, @] ■ 



The study of oxide heterostructures is rapidly devel- 
oping into one of the most promising areas of research 
in strongly correlated electronic systems. The current 
excitement in this field was in part triggered by the re- 
cent discovery of conducting interfaces, with a substantial 
high carrier mobility, between two insulating perovskites 
■ These results were obtained by growing abrupt lay- 
ers of the insulators LaTi(3-|-)03 and SrTi(4-f)03. When 
the spatial distribution of the extra electron was observed 
with an atomic-scale electron beam, it was found to corre- 
spond to a metallic state at the interface [l[ . Theoretical 
investigations of these systems 0, 0, Hi, using Hartree- 
Fock and DMFT techniques, concluded that the leakage 
of charge from one layer to the other explains the results. 
Thus, charge transfer between materials in heterostruc- 
tures can be used to stabilize interface states that other- 
wise would only be obtained via the chemical doping of 
the parent oxide compound. This last procedure has the 
concomitant effect of introducing quenched disorder into 
the sample. Then, for the preparation of doped materials 
without the extra complication of disorder, for the real- 
ization of two-dimensional versions of perovskites, and for 
the exploration of "oxide electronics" devices and func- 
tionalities, the experimental and theoretical analysis of 
oxide heterostructures is a promising and vast field, and 



Before the recent developments mentioned above, fer- 
romagnetic/superconducting (FM/SC) heterostructures 
had already received considerable attention. The analy- 
sis of spin injection from a ferromagnetic material to an 
electrode is of much importance for further progress in 
the area of spintronics. Other properties of interest in 
FM/SC heterostructures include: (i) spin-mixing effects 
that can induce spin-triplet pairing correlations at the in- 
terface |9|] ; (ii) Josephson couplings between two singlet- 
superconducting layers separated by a half-metallic fer- 
romagnet [lo| ; (iii) oscillations in the singlet pairing am- 
plitude detected in SC/FM/SC geometries pT]; and (iv) 
ma gnet ic exchange coupling in FM/SC/FM geometries 
13, 3] ■ Since transition metal oxides have similar lattice 
spacings, the combination of colossal magnetoresistance 



(CMR) manganites [s, 14. Il5l| and high- Tc cuprates has 
been specially investigated [l6l , [itI [18| . Experimental re- 
sults suggest a strong FM/SC interplay resulting in the 
injection of spin polarized carriers into the SC layers 
Of particular importance for the purposes of our investi- 
gations, the transfer of charge from the Lai_a;Caj;Mn03 
(LCMO) ferromagnet to the YBazCuaOy (YBCO) su- 
perconductor has been observed in recent experiments 
using EELS techniques [3]. In this case, the addition of 
electrons, namely the suppression of holes, in the YBCO 
component led to the concomitant suppression of super- 
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conductivity and a transition to an insulating state. Sim- 
ilar effects are also believed to occur in grain boundaries 
of superconductors In other LCMO/YBCO super- 
lattice setups, the superconducting critical temperature 
was also found to be suppressed [21|] . This may be caused 
by the influence of spin-polarized carriers moving from 
LCMO to YBCO, breaking singlet Cooper pairs, or by 
the reduction of the number of hole carriers [19| , or both. 

It is the main purpose of the investigation reported 
here to analyze whether the previously observed elec- 
tronic charge transfer at interfaces from manganites, such 
as LCMO or Lai-^^Sr^^MnOa (LSMO), to hole-doped 
cuprates, such as YBCO, can be extended to the case 
where the cuprate is in an undoped state. If this still 
occurs, namely if estimations of work functions lead us 
to believe that the Fermi level of the manganite remains 
above the first empty states of the undoped cuprate, then 
this Cu oxide will receive an excess of electrons, leading to 
the electron doping of an antiferromagnetic (AF) state at 



heterostructure interfaces 221 . This mixing of cuprates 



and manganites could lead either to (i) a metallic spin- 
polarized state at the interface, if the manganite is fer- 
romagnetic, or to (ii) an electron-doped superconductor 
at low enough temperatures, if the spin polarization is 
not strong enough to destroy Cooper pairs or if the man- 
ganite used is not ferromagnetic. This mechanism would 
allow for the electronic doping of a parent cuprate com- 
pound without the complication of the Coulombic and 
structural disorder introduced by chemical doping. The 
fact that this is electron doping, as opposed to the al- 
ready investigated interfacial hole doping, is particularly 
interesting since a variety of material issues has prevented 
the electron-doped branch of high- Tc superconductors to 
develop as vastly as those that are hole doped. 

Once again, note that if the manganite is ferromag- 
netic, then the charge transfer to the undoped cuprate 
would occur via polarized electrons, which suppress spin 
singlet superconductivity. Thus, there are two competing 
tendencies: adding electrons favors superconductivity, 
but its spin polarization suppresses it. The outcome is 
difficult to predict. Recent investigations using realistic 
microscopic models focusing on a manganite/insulator 
interface predict that a complex pattern of phase sep- 
aration between FM and AF regions could occur [23|. 
Thus, even if superconductivity is not induced, the re- 
sulting state would be interesting, due to the inter- 
play between ferro and antiferromagnetic tendencies. 
More importantly, note that the manganite involved 
in the heterostructure does not need to be ferromag- 
netic: low-bandwidth Mn oxides present a wide variety 
of spin/charge/ orbital order states, other than ferromag- 
netism 

0, 

Il5| . Even large-bandwidth manganites have 
non-FM states at small and large enough hole concen- 
tration. Thus, inducing a SC state is a real possibility 
with the proper choice of the manganite partner to the 
cuprate. 



A qualitative representation of the transfer of charge is 
shown in FigJ^for the cases of a doped cuprate, such as 
discussed in [l9|, and for an undoped cuprate, the focus 
of our investigation reported below. In case (a), SC is de- 
stroyed, while in the proposed case (b), SC is induced (if 
the manganite is not ferromagnetic). Our main result is 
that the novel scenario (b) appears possible. Our conclu- 
sion is based on (i) a detailed analysis of the work func- 
tions of the materials of relevance for these investigations, 
mainly involving the study of photoemission and diffu- 
sion voltage published data, as well as (ii) rough calcula- 
tions at the mean-field and dynamical mean-field levels, 
showing that there is no fundamental problem prevent- 
ing superconductivity from happening at the interface. 
However, it must be clearly expressed that our contri- 
bution should be considered as only the first theoretical 
steps toward realizing electron doping at cuprate inter- 
faces, and considerable more work remains to be carried 
out. For instance, issues such as lattice reconstructions 
and polarity of the involved interfaces have not been an- 
alyzed here (only electronic reconstructions at idealized 
perfect interfaces were investigated). Vacancies and the 
most optimal location of oxygens at the interface have 
not been investigated as well. However, even with these 
caveats we believe the conclusions described below are 
sufficiently interesting that they deserve experimental ef- 
forts to attempt their realization in real heterostructures. 
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FIG. 1: Schematic representation of the systems studied here, 
(a) is a heterostructure involving a ferromagnetic (FM) half- 
metal system, such as LCMO, and a superconductor (SC), 
such as YBCO. In this case, an antiferromagnetic (AF) insu- 
lating state is induced at the interface, (b) is a heterostructure 
involving a non-FM manganite and an undoped AF cuprate. 
Here, the transfer of charge may lead to a superconducting 
state at the interface. Arrows indicate the direction of flow of 
electrons, which are accumulated in the gray region. 



The organization of the manuscript is as follows: first, 
we carry out a detailed analysis of available experimental 
information to judge if the Fermi level of manganites is 
above the lowest-energy empty state of several undoped 
cuprates. This is described in detail, because it provides a 
systematic procedure to address other oxide heterostruc- 
tures in the future. Second, a model calculation in the 
mean- field approximation is carried out, both for static 
and dynamical cases. No fundamental problem is found 
with the proposal of having superconductivity induced at 
the manganite/cuprate interface. The paper concludes 
with a discussion and summary. 
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ESTIMATION OF OXIDES CHEMICAL 
POTENTIAL DIFFERENCES USING 
PHOTOEMISSION AND DIFFUSION VOLTAGE 
EXPERIMENTS 

Overview 

For the success of the emergent field of oxide het- 
erostructures it is crucial to properly determine the rela- 
tive work functions of the materials involved, since these 
work functions control the curvature of the valence and 
conduction bands (VB and CB) of the constituent ma- 
terials, and ultimately the carrier concentration at the 
interfaces. Work functions of conventional metals and 
semiconductors have been studied for decades establish- 
ing the fundamental background for current electronics. 
Thus, for the next-generation electronic devices utiliz- 
ing the complex properties of correlated-electron sys- 
tems, such as high-Tc cuprates and CMR manganites, 
determining the relative work functions of a variety of 
transition- metal oxides is equally important. 

Photoemission spectroscopy (PES) is an important 
technique in this context. PES has provided fundamental 
information to uncover the properties of complex oxides. 
However, although PES techniques can be of considerable 
help for oxide heterostructures, only a limited number of 
experiments have been reported specifically addressing 
the work function of transition-metal oxides 24 1 . 

By measuring the diffusion voltage Vd (or built-in po- 
tential) of a junction between two materials of interest, 
one can also extract the chemical potential differences 
between them. This diffusion voltage is the potential 
barrier at the interface after the rearrangement of charge 
occurs. If there are no extra contributions to Vd, such as 
those caused by interface polarities, impurities, or lattice 
reconstructions, then Vd is equivalent to the work func- 
tion difference between the two constituents. If the work 
function of one of the materials is known, the work func- 
tion of the other can be estimated. This provides another 
procedure to study the band alignment of oxides. 

The purpose of this section is to provide a rough es- 
timation of the band diagram of perovskite transition- 
metal oxides (cuprates and manganites) using the ex- 
perimental data currently available. This is achieved by 
combining information from chemical potential shifts ob- 
tained using PES with diffusion voltage measurements 
on heterostructures. 



Parent compounds of High-Tc cuprates 

Since the discovery of high-Tc superconductivity in 
the cuprates, considerable experimental information has 
been accumulated about these compounds. In particular, 
the chemical-potential jump from the hole-doped to the 
electron-doped materials has been one of the important 



topics since it is closely related to the Mott /Hubbard (or 
charge transfer) gap. Here, let us first consider La2Cu04 
(LCO) and Nd2Cu04 (NCO) as the parent compounds 
of hole-doped and electron-doped cuprates, respectively. 
It should be mentioned that these compounds have dif- 
ferent structures: LCO has the layered perovskite or 
ri = 1 Ruddlesden-Popper structure, while NCO presents 
the so-called T'-structure without apical oxygen. Even 
though the electronic properties of the Cu02 planes are 
quite similar for LCO and NCO, the different overall 
structures suggest that small differences may exist in 
their chemical potentials, as discussed below. 

From the optical measurements, it has been suggested 
that the chemical potential of LCO is located 0.4 eV 
above the top of the valence band 1251 . Similar results 
were obtained with PES techniques [26| . Polaronic effects 
are believed to be the cause of this somewhat "mysteri- 
ous" shift between the top of the valence band and the ac- 
tual location of the chemical potential in LCO ^T^. Sim- 
ilarly, resonant-photoemission studies on NCO revealed 
that the chemical potential of NCO is located 0.7 eV 
above the valence band 2^ 2^. If the positions of the 



valence bands are identical between LCO and NCO, then 
their chemical potential jump becomes 0.7 — 0.4 = 0.3 eV 
(30| . Although an accurate determination of the individ- 
ual chemical potentials of these materials remains to be 
done, estimations based on the shifts of the O Is and Cu 
2p core levels indicate that the chemical potential jump 
is at most 0.5 eV [IH. Then, the positions of the valence 
bands can be assumed to be the same, although with a 
^0.2 eV uncertainty. 

To complete the band diagram, information about the 
location of the unoccupied conduction, or upper Hubbard 
band, is needed. From optical spectroscopy, it was first 
suggested that the Mott gaps of LCO and NCO are about 
2 eV and 1.5 eV, respectively However, numerical 
studies of the Hubbard model on finite-size clusters using 
appropriate hopping amplitudes for the cuprates revealed 
that the Mott gap in Cu oxides was in fact an indirect one 
[3^ : the top of the lower- Hubbard band is at momenta 
(±7r/2, ±7r/2) while the bottom of the upper-Hubbard 
band lies at (7r,0) and (0, tt). Thus, the Mott gap is in- 
direct and smaller than the optical gap. The indirect 
nature of this gap has been confirmed by a recent reso- 
nant x-ray study [l^l . Considering the ratio between the 
magnitudes of direct gap and indirect gap from the the- 
oretical considerations and the optical gap obtained by 
the experiments, the separation between the VB and CB 
of LCO and NCO is estimated to be about 1.5 eV and 
1 eV, respectively j35| . Combining all this information 
together, the schematic band diagram for these materials 
can be constructed and it is shown in Fig. [2] 

Next, let us address another series of high-Tc cuprates: 
YBCO. The diffusion voltage of the heterostructure be- 
tween YBCO (oxygen contents unclear) and 0.05wt% Nb- 
doped SrTiOa (STO) (Nbo.o5-STO) has been measured 
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FIG. 2: Schematic band diagrams of NCO and LCO ^ based 
on the chemical potential jump between NCO and LCO found 
in a resonant-photoemission study [2^ , revised by considering 
the recent discovery of an indirect Mott gap in the cuprates 
[3^ . Chemical potentials (dashed lines) of NCO and LCO 
are located about 0.7 eV and 0.4 eV above the valence bands, 
respectively. Note that the top of the valence bands of NCO 
and LCO do not necessarily match, but estimations discussed 
in the text locate them very close to one another. There 
remains an ambiguity in the chemical potential jump between 
these two materials, but it is believed to be at most 0.5 eV 
am. For more details see text. 



[36|. The chemical potential of YBCO was found to be 
1.5 eV below that of Nbo.o5-STO. 

To establish a connection with the single-layer 
cuprates, note that similar experiments on the diffusion 
voltage of the single-layer parent compound Sm2Cu04 
(SCO) and 0.01wt% Nb-doped SrTiOWNbo.oi-STO) 



have been very recently performed [37|. SCO is 
an isostructural material of NCO, with the same T'- 
structure and formal valences, therefore NCO is expected 
to have a chemical potential very similar to that of SCO. 
The recent experiments of Nakamura et al. revealed that 
the chemical potential of SCO is 1.3 eV below that of 
Nbo.oi-STO 37|. Although there is a slight difference in 
the doping concentration of Nb used for the two exper- 
iments (0.05wt% Nb for YBCO and 0.01wt% for SCO), 
by reducing the Nb concentration to 0.01wt% the diffu- 
sion voltage is expected to increase only by 0.1 eV or 
less because of band-gap narrowing effects 38|. Thus, by 
combining the diffusion voltages on SCO- and YBCO- 
based heterostructures involving a Nb-doped STO sub- 
strate, the chemical potential of SCO is estimated to be 
~ 1.5 + 0.1 - 1.3 = 0.3 eV above that of YBCO. 

Furthermore, assuming that the behavior of SCO and 
NCO are identical, based on the similarity of their struc- 
tures, then the chemical potential of LCO can be esti- 
mated by considering the chemical potential jump be- 
tween NCO and LCO. Figure [3] summarizes the band 
diagram for the various high- Tc cuprates studied here. 
From these results, it is deduced that the parent com- 
pounds of hole-doped cuprates have fairly similar chem- 
ical potentials despite their different crystal structures. 



Doping dependence of the chemical potential in 
High- Tc cuprates 

The doping dependence of the chemical potential is im- 
portant when using doped compounds for the oxide het- 
erostructures. The chemical potential shifts from their 
parent compounds of various high- Tc cuprates have been 
intensively studied 2^, 41, 42, ij]. For the benefit of the 
readers, the main results are summarized in Fig.|4](a) for 
hole-doped compounds, and Fig.[4](b) for electron-doped 
compounds Q. 

There are still considerable discussions on the 
electronic properties of the underdoped region of 
La2-2:Sr2:Cu04 (LSCO), involving concepts such as 
stripe formation, phase separation, and spin glass. Those 
may explain the flat chemical potential shifts of LSCO at 
< X < 0.15. However, focusing on the overdoped region 
where the mixed-phase complexity is reduced, then the 
chemical potential shifts behave linearly with hole doping 
with a slope of about 0.2 eV per 0.1 carrier concentration, 
which is close to that of electron-doped Nd2-xCea;Cu04 
(NCCO). This may indicate that the band structures 
of the Cu02 planes in those compounds are similar, as 
widely believed. On the other hand, the A/i vs. car- 
rier concentration slope for Bi2Sr2CaCu208+y (Bi2212) 
is larger than for the single-layer compounds, perhaps 
due to its bilayered nature. 



Connecting the High-Tc cuprates with the CMR 
manganites 

As already discussed, both from the fundamental 
physics as well as the device engineering perspective, het- 
erostructures involving high-Tj, cuprates and CMR man- 
ganites are important. Here, we consider the cubic per- 
ovskite manganites first, and then turn to the double- 
layered manganites. 

It should be mentioned that the direct measurement 
of the work function of cubic manganites has been per- 
formed by several groups using PES [H, |4^ and the 
Kelvin method [i^l , and results are approximately consis- 
tent with one another indicating an intrinsic work func- 
tion ~ 4.8 eV for Lai_a;Sra;Mn03 with x = 0.3,0.4. 
These experimental data are important in considering 
interfaces with the cuprates. However, as the individ- 
ual work functions of many cuprates are not yet avail- 
able but only their differences, then we have to con- 
sult other experiments in order to establish a connec- 
tion between manganites and cuprates. Here, we con- 
sider PES experiments that measured the chemical po- 
tential shift in doped manganites [3], and those that 
provided the diffusion-voltage on Lan.sSrn.2MnO.'^ (20% 
hole-doped LSMO)/Nbo.o5-STO heterostructures 

According to the experiment by Matsuno et al. 
(iot . the chemical potential of Lai_2;Sr2;Mn03 (LSMO) 
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FIG. 3: Schematic band diagrams of LCO, SCO(NCO), Nbo.oi-STO, Nbo.os-STO, YBCO, and LSMO based on diffusion 
voltage measurements [36l. IstL Is^. l39j and photoemission spectroscopy [29l. l3ll. l40j] . Tops of valence bands (VB) and bottoms 
of conduction bands (CB) are indicated by solid lines, while chemical potentials are indicated by dashed lines. 




Hole concentration p/Cu Doping Level/Cu X, Sr concentration 



FIG. 4: Chemical-potential shifts of various transition-metal oxides, reproduced from photoemission experiments, (a) 
are the resuhs for hole-doped cuprates. Na-CCOC: Ca2-a;Na;i,Cu02Cl2 [111, LSCO: La2-:rSr:rCu04 H, and Bi2212: 
Bi2Sr2Cai_2,Ra;Cu208-t-!/ (R=Pr, Er) The absolute values at the undoped origin of each series are assumed to be the 

same. Reproduced from [43]. (b) Chemical potential shift of electron-doped cuprates NCCO: Nd2-2^Cea;Cu2 04. Reproduced 
from 44]. (c) Chemical-potential shift of LSMO. These shifts are measured at liquid-nitrogen temperature except for LaMnOa. 
Reproduced from [iol ]. 



changes rather monotonically from the parent compound 
LaMnOa to a value A/z=— 0.5 eV for Lao.4Sro.6Mn03 
(60% hole-doped LSMO) (see Fig. 11(c)). Furthermore, 
diffusion voltage measurements on a Lao.8Sro.2Mn03 
(20% hole-doped LSMO)/Nbn 05-STO heterostructure, 
done by Muramatsu et al. [39|, reported a value Vd = 
0.65 eV. Combining these experiments, the bands of 
LSMO can be aligned relative to those of the cuprates 
obtained in the previous section. Figure [3] summarizes 
the band alignment of cuprates and manganites, which is 
our main result of this section. In particular, the chemical 
potential shifts between the various cuprates of potential 
relevance and the 20% doped LSMO are summarized as 



follows: (i) for SCO A/i = fJ-LSMO - fJ'SCO ^ 0.55 eV, (ii) 
for LCO A// — fiLSMO — f^LCO ^ 0.85-1.05 eV, and (iii) 
for YBCO A^ = ^mlsmo — fiYBCO ~ 0.85 eV. In partic- 
ular, it is natural to expect that electrons be transfered 
from LSMO to the valence band of hole-doped YBCO at 
an interface between the two materials. This has already 
been observed experimentally [l^ . thus it is reassuring 
that the simple procedure followed here is consistent with 
available experimental information [48l]. 

Consider now another class of Mn oxides, the double- 
layered manganites with the so-called n = 2 Ruddlesden- 
Popper structure. Direct measurements of work function 
using photoemission have actually been performed on 



6 



Lai.2Sri.8Mn207- It has been reported that the work 
function of Lai.2Sri.8Mn207 increases with decreasing 
temperature across the Curie temperature (125 K) from 
~ 3.5 eV at 180 K to ~ 3.56 eV at 60 K [H, while 
the simple double-exchange model predicts the opposite 
trend. This fact may indicate the importance of sev- 
eral ingredients neglected in the double-exchange model, 
such as electron-electron interaction, electron-lattice in- 
teraction, and orbital degeneracy. The important point 
worth emphasizing here is that the work function of the 
double-layered manganite is more than 1 eV smaller than 
that of cubic LSMO (whose work function is ^ 4.8 eV) 



461 . 147l | , which means that the chemical potential of 
the double-layer manganite is more than 1 eV higher 
than the cubic perovskite LSMO. Thus, at least naively, 
it should be possible to inject electrons to the undoped 
cuprates using the double-layer manganites. However, 
there are issues that remain to be clarified for this conclu- 
sion to be valid, particularly the role of the surface condi- 
tion: for double-layered manganites, the cleaved surface 
used in PES experiments is always A0^~ (A: A-site ion, 
such as La^+ and Sr^+), thus positively charged. The low 
work function of bilayered manganites could be caused 
by this charged surface. On the other hand, for cubic 
manganites the surface can be either positively charged 
(AO termination) or negatively charged (Mn02 termina- 
tion). For this purpose, photoemission experiments on 
cubic manganites with controlled surface conditions are 
highly desirable. 

Summarizing, in this subsection we discussed the possi- 
ble electron doping of cuprates from manganites at ideal 
interfaces. In the paragraphs above, it was concluded 
that the parent compounds of electron-doped cuprates, 
such as NCO and SCO, are the best candidates for this 
purpose since they have the lowest conduction band (the 
CB in LCO is about 0.5 eV higher than that in NCO, and 
the CB in YBCO is expected to be even higher than those 
of NCO and LCO). To dope electrons into a cuprate par- 
ent compound, the minimal condition is that the chemi- 
cal potential of manganites be higher than the bottom of 
the conduction band (or the chemical potential if they are 
not the same) of the undoped cuprate involved. This con- 
dition seems to be satisfied for cubic manganites LSMO 
in a substantial doping range x that includes FM and 
AF states, when mixed with NCO or SCO Possible 
improvements may be achieved by using double-layered 
manganites, which have about 1 eV higher chemical po- 
tential than cubic manganites. 



SUPERCONDUCTIVITY AT THE 
MANGANITE/CUPRATE INTERFACE: STATIC 
MEAN-FIELD APPROXIMATION 

The simple ideas described in the previous sections on 
the possibility of electron doping of cuprates at interfaces 



need to be checked in more theoretical detail to confirm 
their consistency. For this purpose, first a static mean- 
field study will be here discussed, followed in the next 
section by a more detailed analysis including dynami- 
cal effects. The emphasis will be given to a mangan- 
ite/cuprate interface, but the generation of superconduc- 
tivity by electron doping can occur for any other combi- 
nation where electrons are donated to the cuprates. 

The model used here is defined by the following Hamil- 
tonian on a three-dimensional cubic lattice: 



L/Rj 



(1) 

Here Cj^ is the electron creation operator at site i = 
{i-XTiy^iz) with spin s =t:J.j is the electron density 
at site i, and is the nearest-neighbor hopping, which 
is t on the xy plane and tz in the z direction, (pi is the 
electronic potential (discussed in more detail later) that 
will take into account effects related to the charge redis- 
tribution. ^ is the chemical potential, and Wl {W-r) are 
site potentials for the left (right) side of the system to reg- 
ulate the transfer of charge. Hj contains the interaction 
terms. On the left side of the system (the manganite) 
the interaction is 



<^ = -'^HEEcL(^)a/3^i/5-^i' 



(2) 



a, 13 



which is the standard "double-exchange" term, a = 
(ctj; , (Ty , CTz ) are Pauli matrices, and Si is a classical lo- 
calized spin at site i (|S'i| — 1) representing the t2g spins. 
On the right side of the system, the cuprate, the standard 
repulsive Hubbard interaction (C/>G) is supplemented by 
a nearest-neighbors attraction {V<Q) that favors super- 
conductivity in the d-wave channel: 



"iT"-ii 



V^UiUy 

(ij> 



(3) 



The number operator is riis = c\^Cis, and (ij) indicates a 
pair of nearest-neighbor sites i and j on the xy plane. It 
is well known that this t-U-V model leads to a rich phase 
diagram, that includes superconductivity [131 . Periodic 
(open) boundary conditions are used in the x and y (z) 
directions for a lattice of size L = L^x LyX Lz sites. The 
focus of the results described below is on Lx=Ly=16, but 
other sizes have been checked confirming that size effects 
are not strong in this study. 

Two comments are in order: (1) A fully realistic 
model for manganites should contain two Cg-orbitals and 
a Jahn- Teller electron-phonon coupling. However, it is 
common practice in this context to use just one Cg-orbital 
for simplicity, since several phenomena are common to 
both one and two orbitals [l3]. Moreover, the focus 
here is not on the CMR regime but only on combining a 
standard homogeneous manganite state with a cuprate. 
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Thus, the electron-phonon (e-ph) interaction, crucial to 
understand the CMR effect, is here neglected for sim- 
plicity. In fact, LSMO does not have a large CMR effect, 
thus the e-ph coupling may not be strong in this mate- 
rial. (2) More severe is the approximation carried out on 
the cuprate side since a mean-field approximation with 
an explicitly attractive force is employed. However, this 
approximation will be relaxed below and the Hubbard 
model will be used without an explicit attraction, at the 
price of having to constraint the study to smaller systems 
than those that can be analyzed in the static mean-field 
approximation. 

To treat the many-body interactions described by 
Hj^\ here a simple mean-field approximation is adopted. 
For the first interaction term, the following replacement 
is introduced, 

Uj2^il^il ^^Yl [("iT)*^!! - ("iT>("u)] ■ 

i i 

For the second term, a standard BCS model approxima- 
tion is used: 

V nirij 

(ij) 
(ij) 

where Ay = (cj^Cji). The effect of V is restricted to gen- 
erate in-plane superconducting correlations. To further 
simplify the analysis and to reduce computational efforts, 
two sublattices are assumed in the xy plane since it is 
known that this model has a tendency toward antiferro- 



magnetism at half- filling: e.g., (?T.i|) — (n^^^) and (nf^^) 
for {-ly^+'^y = f (A-sublattice) and ~1 (B-sublattice) , 
respectively. More details can be found in [5T|. 

To properly describe a stable junction made out of ma- 
terials with different electronic densities, it is crucial to 
include long-range Coulomb interactions between elec- 
trons and also with the background of positive charges. 
In this paper, the long-range Coulomb interactions are 
treated within the Hartree approximation by using the 
following Poisson's equation |52| : 



,(B) 



-a [{rii) 



(4) 



where a = e/ea (e, e, and a are the dielectric con- 
stant, electronic charge, and lattice spacing, respec- 
tively), is the background positive charge density 
at site i. In the following, e and a are set to be I. Fur- 
thermore, is kept constant to a number ('^+) on 
the left (right) side of the system. The Poisson's equation 
(|4l) is solved numerically using symmetric discretizations 
in the x and y directions, and a forward discretization 
in the z direction, e.g., (P<j)i/dz^ — 0i+2z — 20i+z + (fii 
(z: unit vector in z direction). The use of the Hartree 
approximation here is not a severe problem: the Poisson 



equation is widely considered to be a reasonable starting 
point to account for charge transfer, and it is much used 
in the study of semiconducting systems. 

To mimic a ferromagnetic half-metallic ground state 
on the left side of the system, the localized spins Si are 
chosen to be ferromagnetically aligned, pointing into the 
z direction. To mimic an A-type AF state, these classi- 
cal spins are chosen with the same orientation in the xy 
planes but opposite between adjacent planes. On the 
cuprate side, the mean-field Hamiltonian Hmf is first 
Fourier transformed to momentum space in the xy plane, 
and the resulting Bogoliubov-De Gennes equation [53| is 
solved numerically by diagonalizing a (ALz x ALz) Hamil- 
tonian matrix for each momentum (kx,ky). The d-wave 
order parameter is Ad(i) = (5(i,x) + S{i, — x) — S{i,y) — 
S{i, -y), where (5(i, j) = (ci|Ci+j| -Ci4Ci+jt)/2, and x and 
y are unit vectors in the x and y directions, respectively. 



FM-SC interface: charge transfer leading to AF 

The description of the numerical results starts with a 
qualitative reproduction of the recent experiments [l9| . 
where transfer of charge from LCMO to YBCO was found 
to destroy superconductivity at the interface. In Fig. [5l 
the case when a FM metal forms an interface with a su- 
perconducting state (as in Fig.[Ha)) is studied. The work 
functions and parameters are chosen such that a trans- 
fer of charge from the left to the right takes place, as 
shown in Fig. Elja) where an accumulation of charge at 
the interface is found. The band-bending picture, which 
appears to be true not only for doped band insulators 
but also for Hubbard insulators 5j| , suggests that there 
will be a finite region where the density will be n=l, 
with the chemical potential moving across the Hubbard 
gap. In Fig. E^b) the local magnetization is shown. On 
the left, it is FM as expected in a one-orbital mangan- 
ite model away from n=0 or I. On most of the right, 
the magnetization is zero as in a regular superconductor. 
However, in the interface region the magnetization devel- 
ops a staggered character suggesting the presence of AF 
order. Thus, as expected from the mean-field formalism 
the density n = 1 is associated with antiferromagnetism. 
Figure[5jc) shows the suppression of the superconducting 
order parameter at the interface. Overall, these simple 
results reproduce properly the experiments where it was 
observed that transfer of charge from a manganite LCMO 
to YBCO led to the suppression of superconductivity at 
the interface Il9l. 



Manganite-Cuprate interface, leading to 
electron-doped SC 

After having crudely reproduced the main qualitative 
aspects of previous experiments, the same formalism can 
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FIG. 5: (Color online) Transfer of charge from a FM half- 
metal to a SC state, inducing an AF interface in the pro- 
cess, to reproduce recent experiments [l^ . Shown are the (a) 
electronic density ni, (b) magnetization mz{i) — nij — nij, 
and (c) d-wave superconducting order parameter Ad(i), as 
a function of layer position for two different sublattices 
[(— 1)'^^*" — the xy plane (denoted by circles and 

crosses). Here i = {ix,iy,iz) is a lattice site position. The 
model parameters used are Jh — St, tz = t, Wh = Wt, 
and n+ = 0.7 for the left side of the system, and U — 4t, 
V = — 3t, tz = O.lt (to simulate weak hopping between the 
Cu oxide layers), Wn = 0, and = 0.7 for the right side of 
the system, a = 1 is set for the whole system which has a size 
i = 16 X 16 X 24. The interface is located at = 12.5. The 
localized spins in the left side of the system are fixed to be 
ferromagnetic. The temperature considered here is very low, 
r=t/400. 



transfer of charge. Finally, in (c) the superconducting 
order parameter is shown to become nonzero at the in- 
terface. Thus, as anticipated from the introduction, an 
electron-doped superconductor is predicted to occur in 
the manganite/cuprate system described here, within a 
simple mean- field approximation, for a pair of oxides with 
the proper location of chemical potentials. 

If the A-type AF manganite state for the left side is 
replaced by a FM state, as before, then a similar transfer 
of charge occurs and within mean-field a SC state is also 
generated. However, it is known that the proximity-effect 
influence of ferromagnetism is detrimental to singlet su- 
perconductivity, thus for experimental realizations of this 
scenario, a non-FM state appears more suitable. 

Note: if the transfer of charge would occur in the other 
direction, namely from the n—1 AF to a doped FM ma- 
terial, then the cuprate interface would have less charge 
and a narrow layer with the properties of a hole — doped 
cuprate could be expected. Thus, the simple motive of 
this effort works both ways at the mean-field level. How- 
ever, the analysis of work functions from the experimen- 
tal viewpoint, already discussed, suggests that in order to 
generate interfacial superconductivity involving a man- 
ganite, electron doping is the most likely outcome since 
the chemical potential of LSMO is above that of LCO, 
NCO, SCO, and YBCO. 



SC at BI-AF and Metal-AF interfaces 



be used to formulate predictions for other systems. Of 
particular interest is the combination of an A-type AF 
state (as it occurs for instance in undoped or highly 
doped LSMO and also in bilayer compounds fllflB^) and 
an undoped cuprate. If the transfer of charge occurs in 
the same directions as before, as already suggested by the 
discussion on experimental work functions, then it would 
be expected that at the interface a density larger than 
n=l would be produced (see the sketch in Fig. [Ifb)). If 
this doping is as large as 5-10%, then an electron doped 
superconductor could be created in a real system. 

The actual calculations are conceptually simple and 
the results are shown in Fig.[Sl In (a), the density is pre- 
sented: once again a robust region in parameter space is 
identified where the transfer of charge occurs from the 
AF manganite to the AF cuprate, leading in this case to 
n larger than 1 at the interface. In (b), the local mag- 
netization is shown. In this case, the G-type AF order 
is seen on most of the right (cuprate) region, but this 
magnetic order disappears at the interface due to the 



The generation of a superconductor via transfer of 
charge from another compound is certainly not restricted 
to occur only when LSMO is involved, but it should hap- 
pen under far more general circumstances. For instance, 
the case of a band-insulator (BI) forming an interface 
with an AF was also studied. Figure [Tl^a) shows the den- 
sity profile, indicating that the parameters of the calcula- 
tion are such that the transfer of charge this time occurs 
from the AF to the BI, inducing a region in the cuprate 
with hole doping, and concomitant superconductivity (as 
shown in (b)). Certainly, electron-doped SC can be in- 
duced as well by adjusting the work functions. 

If instead of a BI, a standard metal is used (modeled 
in our study by merely using Hubbard U^O in a tight- 
binding model), then for appropriate work functions also 
a charge transfer is to be expected (see (c) and (d)). 
Thus, the most important aspect of the problem is to 
identify materials with the proper relative location of 
work functions such that the charge transfer occurs in 
the proper direction, and also such that there is a good 
matching between lattice spacings to avoid generating 
extra complications for the transfer of charge to occur. 
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FIG. 6: (Color online) Transfer of charge from an A-type AF 
state (as it occurs in some doped manganites p^) to an AF 
insulator (LCO, SCO, NCO, or YBCO), inducing an electron- 
doped SC state at the interface. The actual model parameters 
used here are Jh = 8t, tz = t, Wl = 14t, and n\ = 0.7 for 
the left side of the system, and U = Ai, V = — 3t, = Q.lt, 
Wr = 0, and = 1.0 for the right side of the system, a — 1 
is set for the whole system, with L = 16 x 16 x 24 being the 
lattice studied. The interface is located at = 12.5. The 
localized spins in the left side of the systems are fixed to be 
antiferromagnetic in an A-type state, and the temperature of 
the study was T=t/400. 
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FIG. 7; (Color online) Generation of a SC state at the inter- 
face between a band insulator (BI) and an AF insulator and 
also between a standard metal and an AF. As in other figures, 
shown are the electron density iii [(a) and (c)], and staggered 
magnetization {—ly'ruzii) and d-wave superconducting order 
parameter Ad(i) [(b) and (d)] as a function of layer position 
iz- The model parameters used in (a) and (b) are Jh ~ (to 
avoid ferromagnetism), = t, Wl = 0, and n+ — 0.0 (to get 
a crude band insulator) for the left side of the system, and 
U ^ 4:t, V = -St, tz = O.lt, Wr = 0, and = 1.0 for the 
right side of the system. The same parameters are used in (c) 
and (d) except nljl = 0.4 for the left side of the system, a — 1 
is set for the whole system with L = 16 x 16 x 24 being the 
lattice studied, and the interface is located at iz — 12.5. 



To place our findings on a firmer footing, we describe 
here the results obtained from cluster dynamical mean 
field theory (CDMFT) [s^ generalized to a layered ge- 
ometry for the case of a ferromagnetic - Mott insulator 
interface. The previous section showed that there are 
many similarities in the use of FM or AF manganites in 
the heterostructure, since their main role is the donation 
of carriers to the cuprate. Since technically the case of 
FM is simpler, in this section we focus on a FM man- 
ganite. CDMFT is a powerful technique which maps a 
full many-body problem onto that of a cluster embed- 
ded in a self-consistent medium. The model Hamilto- 
nian considered here is identical to Eqs. 1 and 2, with 
the addition of the local Coulomb interaction represented 
by the Hubbard U on the Mott insulating side of the 
interface (and without the explicitly attractive nearest- 
neighbors attraction) . Long-range Coulomb interactions 
are treated again using the Poisson equation as in Eq. 4. 

CDMFT can be generalized to treat the effects of a 
layered geometry in the same manner as DMFT was 
adapted to treat a film geometry 0, [s^ . Each layer is 



mapped onto an independent cluster impurity embed- 
ded in a medium. The solution for the self-energy of 
each layer is then used via a self-consistency condition 
to obtain the local Green's function of the entire lattice. 
Hopping between the layers is taken into account only at 
the level of the self-consistency condition. To make the 
problem computationally feasible, we fix the number of 
layers to Lfm = 5 and Lmi = 5 for the ferromagnet and 
the Mott insulator, respectively. Further, each layer is 
mapped to a 4-site cluster embedded in a bath of 8 sites. 
SC correlations are included in the bath around the MI 
layers to allow for the SC instability and a paramag- 
netic solution is imposed. Individual layers are oriented 
in the xy-plane and stacked in the z-direction. The al- 
gorithm starts by making an initial guess for the local 
chemical potential and the bath parameters across the 
layered structure. Solving the cluster impurity model 
in each layer using the Lanczos method we obtain the 
cluster Green's function, cluster density, and self-energy. 
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The density in each layer is used to obtain the new lo- 
cal potential (j)z using the Poisson equation (Eq. 4). The 
self-energy in each layer is used to obtain a local Green's 
function for the lattice using the equation below: 

Giociz^u;)^ [ '^^^^^G{z,z\K.,,Ky,uj) (5) 



where 

G{z,z' ,K^,Ky,uj) = 

[CJ + - t{Kx:,Ky,Z, z') ~ E(z, z', Kx, K, 



(6) 



Here, t^K^^ Ky, z, z') denotes the Fourier transform of the 
hopping matrix (both intralayer and inter layer), and 
Ky denote the superlattice momenta within each layer, 
and z denotes the layer index, /z^ includes all the terms 
that couple to the density, including the chemical poten- 
tial, </)z, site potential, and the work function on either 
side of the interface. The Green's function Gioc{z,uj) is 
combined with the Dyson equation to give a new Weiss 
field that is then used to obtain a new set of bath param- 
eters using a conjugate gradient minimization. Conver- 
gence is achieved when the bath parameters for each layer 
and the density profile across the layers do not change 
with further iterations. 
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FIG. 8: Transfer of charge from a FM manganite to a Mott 
insulating cuprate inducing a d-wave SC state at the inter- 
face in a 10-layer heterostructure, obtained using CDMFT 
generalized to a layered geometry, (a) is the electronic den- 
sity profile Tiz; (b) is the magnetic moment in each layer z, 
denote by mz{iz)', and (c) is the d-wave SC order parameter, 
) across the heterostructure. The model 



parameters used are Jh = 9>t, = t, Wl 
= 1.0, = 0, Q = 1, and C//f=8. 
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In Fig. [5]Ja) the density profile is presented. The bulk 
density is at n=0.7 on the left side (Mn oxide) and ri=1.0 



on the right (Cu oxide, U/t=8). As expected from the 
previous discussions, electron doping of the Mott insu- 
lator occurs at the interface due to a transfer of charge 
from the manganite to the cuprate. Figure [8ljb) provides 
the magnetization at each layer. Note that the effects of 
penetration of the FM moment onto the Mott insulator 
side are not considered in this calculation. In a real setup, 
it is likely that a finite polarization will be found on the 
right hand side of the figure. Figure [D^c) shows the d- 
wave SC order parameter, which in this heterostructure 
becomes non-zero in the cuprate's interfacial region and 
decays rapidly into the bulk of the Mott insulator. Thus, 
an electron-doped d-wave superconductor can be realized 
at the interface of a manganite and a cuprate, due to the 
expected transfer of charge between them, even when dy- 
namical effects in the mean-field approximation are taken 
into account. 



DISCUSSION AND SUMMARY 

In this manuscript, the possible charge transfer from a 
manganite to an undoped cuprate was discussed in the 
context of oxide heterostructures. This issue is nontrivial 
since a priori the existence of a robust gap in the un- 
doped Cu oxides would have suggested the lack of avail- 
able states for manganite electrons to pour into the anti- 
ferromagnetic cuprates near the interfaces. The recently 
discovered indirect nature of the Cu-oxides gaps effec- 
tively reduces the magnitude of the Hubbard gap. Tak- 
ing this into consideration, our analysis suggests that the 
Fermi level of manganites, in a robust range of hole dop- 
ing, lies above the chemical potential of several cuprates 
and even above the bottom of the upper Hubbard band 
of SCO and NCO. This opens the possibility of electron 
doping of high-Tc's at interfaces with manganites (and 
other oxides). Under idealized conditions, the doping 
discussed here may induce a superconducting state (if 
the manganite is not ferromagnetic), an exotic metallic 
state containing polarized carriers in an antiferromag- 
netic background (if the manganite is ferromagnetic), or 
still a superconducting state, likely with a triplet compo- 
nent, if the spin polarization of the carriers coming from 
the manganite is only partial. 

An important component of our effort is the introduc- 
tion of a systematic procedure to analyze photoemission 
and diffusion voltage experiments to predict the direc- 
tion of flow of charge at interfaces. A figure with the 
relative Fermi levels and gap locations of several oxides 
was presented. 

If the superconducting state is ever realized in clean 
interfaces, a fundamental issue to address is the unveil- 
ing of the phase diagram of cuprates in the absence of 
quenched disorder. This is true for both electron and 
hole doping. Recent phenomenological calculations 51 1 



involving noninteracting carriers in interaction with the 
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AF and SC order parameters (an extension of the tra- 
ditional Landau-Ginzburg approach) suggest that these 
phases should be separated by either a region of local 
coexistence or a first-order transition. A glassy state, 
as in the widely studied phase diagram of LSCO, was 
reproduced only incorporating quenched disorder. For a 
sketch of these diagrams see Fig. [9l The generation of su- 
perconductivity at interfaces may reveal the true phase 
diagram of clean cuprates. 




A 


(b) 


AF 






SC \ 




1> 




FIG. 9: (a) and (b) are the possible phase diagrams of the 
cuprates in the clean limit [5^, which may be experimen- 
tally realized at the interfaces discussed in this paper. No 
distinction is made between hole or electron doping, x repre- 
sents both, (c) is the well-known phase diagram of chemically 
doped LSCO. According to Ref. [5l|], the glassy state between 
AF and SC phases is caused by quenched disorder. 



It is clear that our calculations can only be considered 
as suggestive, and its main purpose is to induce further 
work in this area. Important simplifications employed in 
the study, including the neglect of lattice reconstructions, 
vacancies, and polarity effects, need to be addressed. Ab 
initio calculations are crucial to clarify these issues. And, 
of course, the experimental realization of the interfaces 
discussed here would provide a definitive answer to the 
proposal of electron doping of cuprates at the interfaces. 
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